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CHAPTER 1 
INTRODUCTION 

In the early days of helicopter development, the blades 
were rigidly attached to the rotor hub. It was observed by 
pioneers like Sikorsky and Juan de la Cierva that their 
helicopters, which were called autogyros, could hover and 
climb vertically. Yet the machine failed to fly forward, and 
in fact it would get out of control and rolled over. Cierva 
was especially puzzled by this phenomenon since he was able 
to fly a model built with flexible spars without problem 111 . 

It was later discovered that the way the rotor blades 
were attached to the shaft hub was responsible for the 
control problem with the autogyro. Cierva’s full-scale 
prototype had the rotor rigidly connected to the shaft with 
wires while the model was made of flexible material. Figure 
1.1 illustrates the velocity distribution of a helicopter 
rotor blade in forward flight. During hover, air speed 
encountered by each blade is the same. As the helicopter 
flies forward, the air speed experienced by an individual 
blade varies depending on the azimuth angle and radial 
position during each revolution. On the advancing side, the 
resultant speed is higher than the retreating side. There is 
a region of reverse flow on the retreating side where the 
forward speed exceeds the speed due to blade rotation. 
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Forward speed 



Figure 1.1 Velocity Distribution during Forward Flight 

When stiffly attached to the shaft, such as the case of 
Cierva’s autogyro, each blade had the same pitch setting, 
thus angle of attack. As the resultant velocity was higher 
on the advancing side, more lift was generated. The 
difference in lift between the two sides resulted in a 
rolling couple. Since this was not compensated for, the 
helicopter toppled. On the other hand, if the blades were 
allowed to move up and down freely, the ones on the 
advancing side would begin to flap upward because of higher 
lift. As it flapped up, the angle of attack was decreased, 
resulting in lower lift. Maximum flapping for a blade hinged 
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at the center of the shaft occurred at the nose of the 
helicopter. Meanwhile, experiencing lower lift, the 
retreating blade started to descend as it rotated towards 
the tail. Its angle of attack was decreased until the 
retreating blade reached the position over the tail, where 
the local velocity attained a mean value. When the blades 
were allowed to flap up and down, the rotor maintained a 
fere and aft tilt position such that the lift distribution 
was balanced and the machine remained stable. 

Flapping is a phenomenon arisen from the need to 
balance the moment produced by aerodynamic, inertial and 
centrifugal forces about the Center of Gravity. The front 
and back, or longitudinal, flapping discussed above is 
caused by asymmetric velocity distribution during forward 
flight. Besides longitudinal flapping, a rotor blade also 
flaps laterally due to coning. In order for the rotor blade 
to achieve an equilibrium position during flapping, its 
angle of attack has to be just sufficient to compensate for 
the speed at each point during a revolution. The amount of 
flapping also depends on the stiffness of the rotor. In 
flight, a rotor with flapping hinges will maintain a 
slightly coned shape with the lift providing the upward 
force and the centrifugal force keeping the blades extended. 

The idea of having hinge in a blade was further 
utilised in the rotor technology. By incorporating a second 
hinge in the plane of the rotor disk, the blade can move 
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freely without causing stress in the rotor root. This hinge 
is known as lead-lag hinge and together with flapping hinge, 
result in the fully articulated rotor used extensively 

nowadays . 

The ability for a helicopter rotor to flap is an 
important factor in both stability and control. The attitude 
of the tip-path-plane determines the thrust vector in space. 
To control a modern helicopter, the pilot changes the pitch 
of the blade cyclically about the feathering axes by tilting 
the swash plate. The angle between the rotor shaft and the 
tip-path-plane is defined as the flapping angle while the 
angle between the shaft and the swash plate is called 
feathering angle. In forward flight at constant speed, the 
orientation of tip-path-plane is fixed in space. This can be 
achieved with more than one trim condition. One of such 
conditions has the tip-path-plane perpendicular to the rotor 
shaft. Flapping angle is, by definition, zero and only 
feathering is used to maintain the thrust vector. On the 
other hand, a different center of gravity position or 
horizontal stabilizer setting can result in the same trim 
condition with the swash plate being perpendicular to the 
shaft. In this case, there is no feathering and flapping is 
the sole means of control. Typically, both flapping and 
feathering are used to realize the required tip-path-plane 

orientation . 

Flapping can be introduced by either the pilot as a 
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mechanism of controlling » helicopter, or by gust 
disturbance. Oust is undesirable because it has 
deteriorating eifects on the handling qualities of a 
helicopter. Handling qualities is perceived by pilot in 
terms of both stability and control. The ability of the 
pilot to hold the attitude in the presence of gust is a 
stability characteristic. Conversely, if it takes a lot of 
efforts for him to change the altitude of the helicopte 
under -indy conditions, the pilot -ill vie- the control as 
sluggish. Therefore, it is critical to minimise the flapping 

caused by gust . 

Norman Ham of MIT has been -orking on the idea of IBC 
(Individual Blade Control)! 21 . T he method employs blade 
mounted sensors to measure the flapping motion parameters 
and applies the processed information to control the root 
pitch of each individual blade. As a result, flapping caused 
by gust disturbance is reduced. The net effect is a tighter 
control over the tip-p.th-pl.ne, resulting in improvement of 

the vehicle's handling qualities. 

This thesis attempts to demonstrate the utilisation of 

rotor flapping in synthesising an IBC system for gust 
alleviation. The objective of this study is to illustrate 
and seek to improve Ham's IBC method. H sensor arrangement 
with two accelerometers mounted on the root and tip of a 
blade is proposed for estimating of flapping states 
feedback control. Equivalent swash plate implementation of 
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IBC is also deliberated. The study concludes by addressing 
the concept of general rotor states feedback, of which the 

IBC method is a special case. 

This thesis is consisted of eight chapters. The 
following is a brief summary of each chapter. The blade 
flapping equation of motion is derived in Chapter 2. Ham’s 
original IBC method and a modified IBC scheme called Model 
Reference ( MR IBC ) are examined in Chapter 3, followed by 
simulation study with ideal measurements and relative 
performances of the two methods in Chapter 4. The practical 
aspects of IBC implementation are presented in Chapter 5. 
Different configuration of sensors and their merits are 
considered. In Chapter 6, the realization of IBC using 
equivalent swash plate instead of direct actuator motion is 
discussed. It is shown in Chapter 7 that IBC is a particular 
case of rotor states feedback. The idea of general rotor 
states feedback is further elaborated here. Finally, major 
conclusions of this thesis are given in Chapter 8. 
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CHAPTER 2 


BLADE FLAPPING EQUATION OF MOTION 
The equation governing the rigid flapping motion P of a 

single articulated blade in a rotating coordinate system may 

f3,4l. 

be approximated by a linear differential equation : 

(3 + a(u, 0)3 + B ( u , 0 ) 3 = C(u,0)6 + W g (u,w,Q) (1) 

where 

P = blade flapping angle 
H = advance ratio 

0 = blade azimuth angle 

6 = collective and cyclic pitch control input 

Wg = Disturbance gust input 
w = frequency of gust 

Q = rotor speed 

The coefficients in the equation of motion are periodic in 
0, and are function of u. The equation is reduced to a 
ordinary linear differential equation with constant 
coefficients in hover, when u = 0. 

The equation of motion can be derived by summing 
moments about the articulated offset hinge, as shown in 
Figure 2.1. Forces acting on an element of mass dm are : 1) 
aerodynamic lift force dL, 2) force due to tangential 
acceleration of flapping (r-e)P x dm and 3) centrifugal 
force due to rotation about the rotor shaft rQ 2 x dm. 
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rQ^ dm 
(r-e)3 


dm 


Figure 2.1 Forces acting on a blade element 


By using harmonic balance technique 153 , the pitch and 
gust inputs can be expressed in assumed harmonic motion 
forms. The root pitch input 6 assumes a harmonic motion of 
steady-state term plus harmonics at rotor rotation frequency 
Q : 


8 = 0q - A ls cosQt - B ls sinQt 

The gust input assumes the form of a sine wave of 
frequency w and first subharmonic (fl - w) and superharmonic 
(Q + w) of the rotation frequency : 

W g = ai sinwt + bj (cos(Q -w)t - cos(fl + w)t) * U (3) 
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Assumptions made in deriving the above blade flapping 
equation of motion are i 

1) The blade is rigid. 

2) Only a single blade is considered, there is no 

interaction with the fuselage. 

3) Flapping is assumed to take place in the pitch 
plane, there is no coupling with the roll or yaw 

axis . 

4) In deriving the equation, any higher order terms 
beyond second of the normalised flapping hinge 
offset € are dropped. 

5) The inflow of air to the rotor disc due to forward 
flight is lumped into the gust term W g . 

6) The rotor speed is constant. 

7) Gravitational force on the blade is neglected, 
only inertial and aerodynamic forces are being 
considered . 
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CHAPTER 3 

IBC ANALYTICAL FEEDBACK CONTROL 
A helicopter, compared to an airplane of similar site, 
provides a smoother ride in the presence of gust. Since the 
rotor blades flap individually to reduce the pitch and roll 
rate, the fuselage and the rest of the helicopter are 
shielded from the disturbance. In contrast, the wings of an 
airplane simply transmit the wing loads, either as 
disturbance or lift, to the fuselage. The result is a 
rougher ride. 

The air space around a helicopter rotor is a very 
complex aerodynamic environment. The frequency response of a 
blade can be divided into two categories according to the 
frequency contents of the excitation. The first one is the 
low frequency range from 0 to 1 Q (rotational speed), which 
includes gust response, blade instability due to flap-lag 
coupling, and ground resonance. The second category is a 
high frequency domain above 1 Q involving blade bending 
stress, vibration, noise and higher harmonic motion. 

The blade flapping mode, which completes a cycle during one 
blade revolution, has a natural frequency of approximately 
1 Q. This mode dominates low frequency effects such as gust 
response and flying qualities. Since blade flapping is a 
means by which the pilot trims and controls the helicopter, 
any unwanted blade flapping due to gust is viewed as an 
external disturbance to be attenuated. 
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Methods of IBC An alysis. 


To best illustrate IBC technique, general methods used 
to analyze systems described by differential equation with 
periodic coefficients will be discussed, followed by the 
particular IBC control algorithm used by Ham. Finally, a 
modified approach termed MRIBC (Model Reference Individual 
Blade Control), is presented for possible improvement of the 
IBC concept. The equation describing the rigid flapping 
motion of a single articulated blade in a rotating 
coordinate system can be approximated by a linear 
differential equation with periodic coefficients . 

(3 + A(p,0)f 3 + B(u,0)3 = C(p,0)9 + W g (p,w,Q) U) 

Many features of analysis and synthesis of linear time- 
invariant system, such as eigenvalues and linear quadratic 
synthesis, can be applied to a linear periodic system, 
according to Floquet theory' 4 ' 61 . The theory states that if 
the solution of a linear periodic system is sampled at 
intervals of one system period, the sampled solution behaves 
as a time-invariant system. As the first step in designing a 
controller for the IBC system, the periodic system is 
transformed into a time-invariant one. Linear quadratic 
synthesis can then be employed to determine a set of control 
gains which optimizes a particular cost function. 
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mm's IBC (Indi v idual B]a<ie Control ! 

Borman Ham of MIT has bean doing research work on IBC 
concept. IBC is an active control in which each blade is 
separately controlled in the rotating tr.me by attaching 
broad-bang el.ctrohydraulic actuators to the swash plate or 
individually to each blade. The control signal is generated 
from sensors mounted on the blades. IBC involves not just 
control of each blade independently, it also provides a 
feedback loop for each blade in the rotating frame. 

In Ham's IBC design, the controller design (blade 
flapping angle, angular rat. and angular acceleration 
feedback gains) was not conducted using multi-variable 
states feedback control methods such as linear quadratic 
analysis. Instead, the feedback gains were selected such 
that the closed-loop blade dynamic behavior was the same 
the open loop's, with the exception of a reduced gust term. 
The derivation of Ham's IBC scheme is given in Appendix A. 
in Chapter 6, it is shown that Ham's IBC soheme is a limited 
case of multi-variable rotor states feedback control system. 
Figure 3.1 shows a block diagram of the open loop system. 

8 S „ is the trimmed collective and cyclic pitch input from 
the swash plate to sustain the forward flight. W, is the 
disturbance gust input. A typical plot of input and output 
for this open loop system is depicted in Figure 3.2. When 
the output angle, rate and acceleration are fed back, as 
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Figure 3.3 Ham’s IBC method 

Periodic controller gains are difficult to implement 
physically. It also places a lot of demands on the actuator. 
A set of constant gains may be chosen for feedback while the 
flapping equation of motion still maintains the same 
periodic coefficients. This configuration is being 
investigated as an alternative to the time-varying gain 
model in a simulation study in Chapter 4. 

Another viable approach to Ham's IBC is apparent when a 
closer look is taken at the equation of motion. The periodic 
coefficients are cosine and sine functions multiplied by 
first and second order terms in y. These terms will not vary 
much when y is close to aero. A constant coefficient model 
of the blade flapping equation may be valid for low advance 
ratio. In reference 2, a Floquet analysis was performed to 
the flapping equation of motion and it was found that the 
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A math model of the plant can be built using technique 
such as parameter identification. The differences between 
the plant and the model states, rather than the plant states 
themselves, are used to formulate the feedback control law. 
The IBC control is a function of the error signal e, and 
therefore less sensitive to DC value of the control input. 
Derivation of the equations for MRIBC method is shown in 

Appendix B. 

Owing to the scope of the present studies, some 
essential subsystems in IBC are not addressed. For instance, 
the actuator dynamics is ignored in the analyses. Actuator 
is basically a low-pass servo system. When it is placed in 
the feedback loop in Figure 3.3, 0 ibc will contain only the 
low frequency portion of the original signal. High frequency 
flapping disturbance will not be affected by the IBC system 
since the high frequency content of the feedback signal is 
filtered. Flapping states (3, (3, (3, used in determining 0 ibc 
are either measured or estimated. In both cases, the signal 
is corrupted due to presence of noise. Low pass filters are 
usually required before this signal can be utilized. This 
will introduce delay (phase shift) in addition to effects on 
the overall stability of the whole IBC system. 
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CHAPTER 4 


SIMULATION ANALYSIS OP IBC WITH IDEAL MEASUREMENTS 

A simulation study was performed to evaluate the IBC 
schemes employing the Model Reference method and the 
original one conceived by Ham. The simulation runs were made 
using a VAX digital computer. A fourth order Runge-Kutta 
method was utilized for numerically solving the differential 
equations of motion. Integration step size was 0.001 second. 
A generic UH-60 Black Hawk rotor blade was chosen for the 
study. Effects due to fuselage interaction or actuator was 
not considered. In the simulation, sensor dynamics was 
ignored and perfect measurements of flapping angle, rate and 
acceleration were presumed available. 

Ham’s IBC Method 

Three approaches were made in analyzing Ham's IBC 
method. The first approach employed the full periodic 
flapping equation of motion with time-varying feedback 
control gains. The second method still used periodic 
coefficients for the equation. The controller gains were, 
however, independent of time. They were computed as average 
value of the time-varying gains. In the third approach, 
constant coefficients were assumed for the flapping motion, 
resulting in constant gains. It is known as simplified Ham's 
IBC method in this study. 

To investigate the effects of IBC on gust alleviation. 
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both the open and closed loop flapping angles were plotted 
on the same graph. Peak to peak values of the two angle 
responses were measured and the percentage of reduction due 
to feedback computed. For example. Figure 4.1 shows! for the 
first approach of periodic coefficients and gains, the 
decrease in flapping angle response when the gain KA takes 
on a value of .35 and 1.2. The dotted curve is open loop 
response (KA = 0) and the solid one represents the response 

when IBC is in effect. 

In Ham's second IBC approach, an effort was made to 
find a set of constant feedback gains which would work 
reasonably well for the time-varying flapping equation of 
motion. It was observed that the system became unstable when 
the gain KA was over 0.55. The results are summarized in 

. - t ^ 

Figure 4.2 for^KA = 0.35 and 0.575. There appears to be a 
phase swift in the angle response between the open loop (KA 
= 0) and KA ^ 0 case. This method is not further pursued in 

the current study. 

Discus sion of Results. 

Sensitivities of flapping angle response to control 
gains are demonstrated in Figures 4.3 and 4.4 for first and 
third approach respectively. Gust effect is seen to be 
reduced by increasing the value of KA. For any given value 
of KA , the first approach with time-varying gains tends to 
achieve a better performance in reducing flapping due to 
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gust disturbance. On the other hand, the simplified IBC is 
easier to implement since a constant set of feedback gains 
is used through out the flight. 

As shown in Figures 4.5 and 4.6, the magnitude of the 
IBC control input, 0^ c > also gets larger with increasing KA. 
Furthermore, even though the magnitude of the total control 
input TTOL , defined as the sum of IBC control and swash 
plate input, remains about the same, its DC value increases 
as KA is changed from 0.35 to 1.2. This increase in DC value 
is equivalent to adding an extra collective pitch input due 
to swash plate motion. An excessive DC value may impair the 
vehicle's flight qualities. The allowable peak to peak value 
of the control input places an upper limit on KA. A 
realistic upper bound of KA is about 1.2. 

From the derivation in Appendix A, it is established 
that, in theory, the gust term can be reduced by a factor of 
1/Kg using Ham's IBC method. For the simplified IBC scheme 
with constant coefficients and gains, Kg - 1 + KA. Thus, a 
theoretical reduction of 50% in gust input can be 
accomplished with KA = 1 for the simplified IBC method. The 
corresponding reduction in flapping response (3 is not 
necessarily the same in magnitude. The simulation results in 
Figure 4.8 show that a reduction of 20% in flapping response 
is achieved when KA = 1 . 

MRIBC (Model Reference IBCl Method 
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In this study, the plant was modelled by differential 
equation with time-varying coefficients. Two model reference 
control schemes, one with time-varying math model and 
feedback gains and the other with time-invariant model and 
gains were simulated. 

Discussion of Results 

A typical time history of ^ibc and 13 for MRIBC with 
time-varying math model is shown in Figure 4.7. Note that 
the feedback signal 0 ibc has a DC value of rero. This is a 
characteristic of MRIBC since £ depends only on the 
differences in the plant and model states. MRIBC control 
activity is minimum with absence of gust disturbance on the 
blade. Between the methods with constant math model and 
time-varying model, the latter exhibits higher percentage of 
gust reduction for the same KA. The results are summarized 
in Figures 4.8 and 4.9. As KA is increased, both models 
become unstable eventually. This happens at KA > 0.75 for 
the time-varying math model . The model with constant 
coefficients will be unstable if KA > 1.2. 

MRIBC with time-varying math model of the plant shows 
performance comparable to that of the original IBC s. In 
fact, for KA = 0.5, MRIBC has a slight edge over IBC. 
However, MRIBC becomes unstable rather quickly. MRIBC with 
time-invariant math model generally has smaller percentage 
of gust reduction than simplified IBC. 
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As previously discussed, MRIBC has an control input 
whose DC value is less sensitive to variation of controller 
gain KA than IBC. This property is illustrated in Figures 

4.10 and 4.11. While Ham’s simplified IBC is shown in Figure 

4.11 to have a $ with DC value of about -0.2 rad., the 
corresponding value is essentially zero in the case of MRIBC 
with time-invariant math model. 
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CHAPTER 5 


ESTIMATION OF BLADE STATE VARIABLES 
For any of the above IBC schemes to work, blade angular 
displacement, angular rate and angular acceleration are 
needed to generate the feedback signal . The information can 
definitely be obtained by direct sensor measurements, as in 
the case of full blade state variables feedback. However, 
there are several disadvantages with measuring all states 
directly. First of all, sensors are expensive and they add 
complexity to the overall system. Secondly, there is always 
noise associated with any measurement. In certain cases, the 
noise present may make an otherwise stable system unstable. 
To counteract the noise, filters are invariably required for 
these sensor data, again adding cost and complexity. 

Figure 5.1 shows a block diagram of a typical closed- 
loop estimator. 



Figure 5.1 Closed-loop estimator 
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The idee of estimator is to make use of algorithm to 
reconstruct all the blade states given only measurement of 
some of them. A plant is represented by a set of first order 

differential equations . 


X = 


Let X = 

A 

X = 


F X + G U 
estimate of X 
F X + G U + K(Y 
F X + G U + K(Y 


Y) 

H X) 


( 1 ) 


( 2 ) 


Let 


= estimation error = X - X 


X 


, A 

X - X 


F X + GU - [FX - GU + K(Y 
(F - K H) X 


HX)] 


This is a homogeneous differential equation. Values of 
the estimator feedback gain vector F can be chosen such that 
the estimator system matrix (F - F H> represent, a stable 
system. Furthermore, F should be selected such that the 

. fVt*r#fore the estimation error, 

estimator system, and theretore rne ^ 

converges to sero as fast as possible. When this happens, X 

will converge to X regardless of the value of X(0>. For the 

selection of the gain vector K, the characteristic equation 

of the system matrix (F - F H) is computed. The coefficients 

of like powers of this equation are then compared to the 

equation formed with the desired estimator root locations in 
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the s-plane. K is uniquely determined if the output variable 
y is a scalar, i.e. there is only one non-zero element m 
the row matrix H. 

Estimator roots are picked to be faster than those of 
the plant so that total response is dominated by the 
response due to the plant. The faster the estimator roots, 
the quicker the error converges. Upper limit to the 
estimator response speed depends on noise rejection 
characteristics and sensitivity to model errors. In the 
simulation, the estimator roots were selected at least 4 
times faster than the plant's. In the following section, the 
performances of two estimators, namely Direct Flap Estimator 
and IBC Kinematic Estimator are evaluated. 

Direct Flap Estimator 

Direct flap estimator uses only the flapping angle 0 to 
estimate 0, 0 and 0* In state space form, the open loop 
flapping equation of motion is written in the form of (1) 
as : 




0 I 


d 


; x = 

• 

-B -A 


0 





u = ®swp + w g 
G = [0, 1] T 

H = [1, 0] 
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In state space form, the estimator for the blade flapping 
equation of motion is given as : 


d 

P 


0 
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X, 
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0 

dt | 



-B 

-A 


a 

+ 

e +w 

_ swp g 



ip - ft 


(3) 


Case 1 - A , Direct Flap Estimator 

The coefficients A, B are being utilised as feedforward 
terms and the estimation error [|3 - fl] as feedback term. 
Ham's simplified model with constant coefficients and 
controller gains was employed to demonstrate the performance 
of the estimator. The estimator roots were set about 4 times 

faster than the fastest plant root. The result is shown in 

, - * i 

Figure 5.2. The dotted and the solid lines represent the 
true and estimated states respectively. The estimated angle 
follows the true one very closely and the rate is doing 
reasonably well. 

Case 1-B. Direc t Flap Estimator 

The plant dynamics change with flight conditions, and 
is therefore time-varying. The flapping equation is a 
function of forward air speed and blade atimuth angle. If 
the estimator roots are to be maintained throughout the 
flight, estimator gains will have to be adjusted 
continuously^ . This can be done by gain scheduling as a 
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function of flight conditions. However, this is not an easy 
task as the combinations of flight conditions are numerous. 
If the plant dynamics are taken out by setting A. B and th 
input matria to aero, a set of constant estimator gain, can 
be utilised through out th. flight. The model becomes known 
as kinematic estimator! 1 '' 1 . The equation below describes 
such an estimator. 


d 

/v 

P 


r 

0 

I 


f 

+ 

K 1 

dt 

/S 

P_ 


0 

0 




_ K 2_ 


By analogy with classical control theory. Equation (4) 
characterises PI (proportional plus integral) control since 
the estimation error is used to update both the angle and 
the rate. The effectiveness of the kinematic estimator is 
depicted in Figure 5.3. It does not track the plant as well 
as the previous estimator in Equation (3) which includes the 

plant dynamics. 

The kinematic estimator described by Equation (4) 
employs (3 to estimate both (3 and (3. Since acceleration (3 is 
also needed for the IBC schemes, the estimator in Equation 
(4) will be extended to provide (3 information. Two cases are 
considered. The first one is essentially two Direct^ap^ 
Case 1-B estimators cascaded together. Figure 5.4 
illustrates this method. The output of the first estimation. 


26 



* * f ' 

3 is used as input to predict 3. It is seen from Pigure 5.5 
that while the angle can be tracked very well, the estimated 
acceleration lags behind the true one. The second case is 
delineated by Equation 5.5. 


d 

dt 


A 

P 

A 

4 

P 

<7 

P 


r i 


A 



o 

H 

O 


P 


K i 






0 0 I 


p 

+ 

K 2 



c 



0 0 0 




_ K 3_ 


IP - PI 


(5) 


The system matrix is 3 x 3. Three estimator roots close 
to the ones used in equation (4) were chosen to evaluate 
this implementation. The performance is similar to case one 
above, as shown in Pigure 5.6. 



Pigure 5.4 Cascaded Direct Plap estimator 
IBC Kinematic Estimator 

Por the Direct Plap estimators described above, 
measurement of flapping angle is required to estimate $, 3 
and 3. This can be done by attaching a strain gauge at the 
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blade root fitting. A better technique to deduce all three 
states needed for IBC is to use two accelerometers mounted 
at the root and tip of a blade™. In this IBC Kinematic 
Estimator method, the two accelerometers provide both (3 and 
p information, from which P can also be estimated. A block 
diagram for this approach is given in Figure 5.7. 



Figure 5.7 IBC Kinematic Estimator block diagram 


An accelerometer positioned along an hinged blade 
experiences inertial force as shown in Figure 5.8. 
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The signal a { (x) is in response to both the flapping 
and first elastic flatwise bending mode, 
af ( x ) is given by 

a f (x) = RQ 2 ( x( 3 /Q 2 + x(3 + n(x)g/Q X + xn(x)g ) ( 6 ) 

where R = rotor radius 

Q = rotor speed 

x = blade spanwise location (x - e)/R 

3 = blade flapping angle 

n(x) = first elastic flatwise bending mode shape 

g = first elastic flatwise bending mode 

displacement 


The first elastic bending mode has a natural frequency 
at about 2 Q. If a low pass filter with bandwidth up to 1 Q 
is used, the flapping mode, which has a natural frequency 
between 0 and 1 Q, will be preserved while the first elastic 
bending mode is attenuated. With such a low pass filter, the 
accelerometer signal becomes : 

a f (x) = RQ 2 ( x(3/Q 2 + xp ) 

This is an algebra equation with two unknowns on the 
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right hand side. With two accelerometers placed at different 


locations on the blade, two such equations are obtained. 

Both flapping angle and acceleration can then be solved for. 
Written in scalar form, Equation ( 4 ) becomes : 


dP * 

— = P + K, (/?-/>) 

dt 1 


dp 

— = k 2 (p - p) 


( 8 ) 


Equation (8) is modified to accommodate the additional 
flapping acceleration information. The IBC Kinematic 
Estimator equation can be expressed as : 

dp * 

“ = P + Kj (/? - /?) 

dt x 

dfl (9) 

— - P * k 2 (/>-/>) 

J a. ** 


• * 

Since both 0 and 0 are available, the task of 

. * 
estimating 0 is not difficult. In fact, the estimated 0 

tracks the true one very closely after a transient duration 

of 0.5 second. 
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CHAPTER 6 

IMPLEMENTATION OF IBC THROUGH SWASH PLATE ACTUATORS 

All the IBC methods discussed so far utilise an 
individual actuator to change the pitch of each blade. This 
is direct control in the rotating frame. IBC is executed 
using these actuators, which rotate with the blade while 
pilot control is achieved using a conventional swash plate 
in non-rotating frame. In this configuration, actuator 
reliability is always a concern, and it may be more 
significant than the simplicity such a system offers. As an 
alternative, by attaching actuators to the swash plate, the 
same degrees of individual blade control can be accomplished 

in a pure non-rotating frame. 

A helicopter with swash plate has three control 
degrees-of-f reedom : collective, longitudinal and lateral 

cyclic. For a helicopter with three blades, the number of 
control degrees-of-f reedom is equal to the number of blades. 
When four blades are present, as in the case of Black Hawk, 
an equivalent swash plate implementation of IBC can also be 
formulated. For a four bladed helicopter, a transformation 
between the rotating coordinate system and the non-rotating 
system is given as : 

Afibc = TRA ^vp 

where A 0^ = [ A 0j , * 9 2 , 0 3 , * e 4 ] T is the individual blade 

pitch control angle vector in the rotating reference frame. 
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and £ fLg W p = [ £ 0 cl , A A ls , £ b 1s , ^ 9} ] T is the ideal swash plate 
displacement vector in the non-rotating frame. is any 
element ofAe^c* It can be expressed in elements of swash 
plate displacement vector as 

£ 0 i = £ e c ] + £ Ais cos *i B ls sin *i (-1) 1 (D 


for i = 1,2, 3, 4 in the case of a four-blade system. 
Thus, the transformation matrix can be expressed as 



1 

C#- 

S#i 

-1 

TR = 

1 

-s#; 

C#i 

1 


1 

-C#T 

-s#i 

-1 


1 

s #i 

-c#i 

1 


( 2 ) 


where # = azimuth angle, C#! = Cos#!, S#! = Sin#! 

An equivalent swash plate motion for a given ®ibc can ke 
derived as 


^ Qgvrp ~ TR * £ 6^bc 
where 


TR 


-1 _ 


1/2 


1/2 

C#i 

S#! 

1/2 


1/2 1/2 1/2 
-S#! -C#! S#! 

C#! -S#! -C#! 

1/2 1/2 1/2 


( 3 ) 


Since the current swash plate mechanism has only three 
control degrees of freedom, does not exist. Only 

[ £ 0 c j ,£A! S A B 1s ] motion can be realired in real application. 
If they are fed back to the rotor system, then 
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( 4 ) 
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1/2 1/2 1/2 
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A i s 
B Is 

= 1/2 
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S0i 

-S 0 ! -C 0 ! S 0 X 
C 0 ! -S 0 ! ”C 0 ] 


©2 

%_ 


The equivalent blade pitch angle becomes 


63 

0 4 
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1 

1 
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S0! 


S0i 

C0! 

-S 0 ! 
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0 cl 

1 


*ls 

-1 
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1 





__ 


When Equation ( 4 ) is substituted into ( 5 ) 




0 1 


0 2 

s 

0 3 


0 4 



3/4 

3/4 

-1/4 

1/4 

1/4 

3/4 

1/4 

-1/4 

-1/4 

1/4 

3/4 

1/4 

1/4 

-1/4 

1/4 

3/4 


0 1 

0 2 

e 3 


e 


4 


( 5 ) 


( 6 ) 


It can be seen from the above relationship that an 

identity matrix between [ ^ 0 i, a 02 ' A0 3' A0 4 and 
f e, , Sj, 63, 64 ] T is not present in general. The two vectors 

in Equation (6) will be identical if ( andAe 3 ) and 

( 9 2 and 0 4 ) have the same magnitude and are always 100 
degree out-of-phase with each other. A simulation study with 
the Black Hawk GENHEL model was conducted and the result is 
depicted in Figure 6.1 The equivalent blade root pitch 

A 

motion due to IBC generated through the swash plate, A 0 ]., is 
reasonably close to the original direct IBC motion , A 0 1 , in 
rotating frame. Given the complexity in mounting the 
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actuators directly on the blades and reliability 
considerations, IBC through swash plate is a practical way 
for implementation of IBC. 



J x 2 B i M » n V* J » j * 


a e 

A ^ 

A 0 


: blade root pitch generated by direct IBC 

blade root pitch generated b> IBC via swash-plate motion 

Figure 6.1 Blade root pitch motion 
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CHAPTER 7 


GENERALIZED ROTOR STATES FEEDBACK CONTROL SYSTEM 
In this chapter, a generalized rotor states feedback 
control system is defined. It is shown that the IBC system 

through swash plate control is a limited case of a general 

[9] 

rotor states feedback control system . 

The general coordinate transformation between blade flapping 
angles H* in rotating frame and the rotor states in non- 

rotating frame is given by : 


fi* = T 

where 

T = 

Hr s [ Hi* &2' 
H-JTR = £ 00' 01c 1 




1 

C0i 

S0^ 

-1 

1 

-S0 : 

C0 X 

1 

1 

-C0j 

-S0 X 

-1 

1 

S0i 

-C0 X 

1 


fl 3 , I3 4 , 1 T 

0 ls , ] T 


( 1 ) 


rotating frame 
non-rotating frame 


(3, = 0 O + (3ic COS0 1 + (3 ls sin0i ♦ (-D 1 

for i * 1, 2, 3, 4 


Generalized pure I B S. 

For Ham’s pure IBC system, individual root pitch angle 
vector AS^c is generated from the blade states C 6 r 0* 0* j 
using sensors mounted on the rotor blades. The control is 
performed in the rotating frame through actuators attached 
to the blade roots. The control law is given as 
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( 2 ) 


^tibc = K X g* + K 2 g* + K 3 % 

Equation (1) is differentiated twice to give 

a* 1 » fim ♦ T (3) 

E* = T am, ♦ 2T en„ + T Sn (4) 

Equation (2) is now written in terms of the non-rotating 
rotor state vector and its derivatives as 

Afiibc = CKi T + K 2 T + Kj T] ♦ 

[2K X T + K 2 T] g#* + 

T] g|(^ ^ 

Equation (5) shows a pure IBC concept through the root pitch 
actuators control is a special case of the rotor states 
feedback control system with time-varying feedback gain 
vector as given in Equation (5). 

Generalized IBC Concent v ia Swash plate Feedback 

For a helicopter with 3 or 4 blades, the IBC control 
can be realized through the swash plate motion. Using 
Equation (1) 


Since 

therefore 


fiibc * T a 






= T' 1 fiibc 


( 6 ) 


A 0 


= T' l [Ki T ♦ K 2 T + K 3 T] firo + 
T‘ 1 [2K 1 T + K 2 T] gjfR + 

T _1 [Ki T] gjR 


( 7 ) 
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This is equivalent IBC via rotor state feedback in the non- 
rotating frame through the swash plate motion. Now, consider 
a simplified case where the control law in Equation (2) is 
given as : 


K]_ = k L I 4 , K 2 = k 2 I 4 , *3 = k 2*4 

Here k^ k 2 and k 3 are scalar feedback gains and I 4 is an 
4x4 matrix. With this simplified control logic. Equation 
(7) becomes : 




Ck x T~*T + k 2 T _1 T + k 3 I 4 ] 0|fR ♦ 
[2k, T' 1 ! + k 2 T‘ L T] + 
t k l ^ fijTR 
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-k,Q 


- kiQ 
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k 3 






k 2 
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•2ki 

0 




0 

2 k,Q 

k 2 
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0 
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0 




k l 
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0 
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k l 

0 
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0 
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k l 
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0 

0 

0 

kl 
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Note only first three control signals (elements) in A0j wp are 
used for gust control . 

For another special case where = K/Q 2 , K 2 = K/Q and 
Kj = K, Equation (7) becomes 

A 9*™ = [K/Q 2 T _1 T + K/Q T _1 T ♦ K] + 

[2K/Q 2 T* 1 T + K/Q T _1 T] fijfR + 


After the matrixes are expanded, Equation (8) becomes 


*S«p 


K 


— _ ; 

1 0 0 0 . 

— “1 
10 0 0 

0 0 1 0 J3^ + K/Q 

0 12 0 

0-1 0 0 j 

0-210 

0 0 0 lj 

0 0 0 1 


£*R 



( 10 ) 


The control shown in Equation (10) is the original IBC 
concept proposed by Ham. Comparing Equation (10) with 
Equation (5), it is cleat that IBC control scheme suggested 
by Ham is a limited case of the rotor states feedback 
control system. 
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nongral Rotor State Feedbac k through the Swash P l at? MpU .211 
All the discussions so far on IBC control do not 
consider the helicopter body. In fact/ motion of the 
fuselage will affect rotor flapping and vice versa. The body 
of the helicopter is isolated from external disturbance 
because of stabilising effect provided by the rotor, which 
acts as a gyroscope. If fuselage body states are available 
for feedback, more versatile gust alleviation control law 
can be formulated. Studies on IBC by J.C. Wang has found 
that a significant reduction on the body vertical 

acceleration caused by gust can be achieved by adding a body 

[9] 

vertical acceleration feedback loop into the IBC system 

If Am can be estimated from sensors mounted on 

the blade, a general gust control law can be developed by 
using rotor and body states feedback as 

0 1 * 

A = K 1 + K 2 £*R + K 3 + K 4 Z* + K 5 
where the fuselage body state vector Xj is 
Xj = [u,q,w,v,p,r] T 


u , v , w : linear velocity 

p,q,r : angular velocity 

The optimal gain matrixes Kj through K5 for a given flight 


39 


condition can be derived for control law synthesis when the 
nine degree-of -freedom rotorcraft state dynamic equations 
become available. 
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CHAPTER 8 


CONCLUSIONS 

In modern day helicopters, flapping provides a means of 
orientating the rotor thrust vector for stability and 
control functions. Simulation studies were performed to get 
some insight into Ham's original IBC method of reducing 
blade flapping due to gust disturbance. Results show that 
for an IBC model with periodic coefficients and controller 
gains, up to 36% in gust alleviation can be achieved with a 
maximum feedback gain KA of 1.2. When a simplified model 
with constant coefficients and gains is used, the reduction 
is decreased to about 22% for the same value of KA. 

It is desirable to have a control system which is less 
sensitive to the DC value of the input control. A MR IBC 
(Model Reference IBC) method is suggested for possible 
improvement of Ham’s original IBC. The sensitivity to 
control input is lessened since differences in flapping 
states between the model and the plant, instead of the 
states, are employed for control law synthesis in MRIBC. The 
performance of this method matches very well with that of 
Ham's in both time-varying and time-invariant cases. The 
Model Reference IBC system, however, becomes unstable at 
KA = .75 and KA = 1.2 respectively for time-varying and 
constant math model . 

For its insensi tivity to flight conditions and 
simplicity of implementation, kinematic estimator is used to 
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reconstruct the states required for IBC feedback qiven only 
measurements of some of the states. Two algorithms for 
estimating both flapping angle, rate and acceleration from a 
single measurement of angle (Direct Flap Estimator) were 
investigated. While satisfactory tracking of angle and rate 
can be attained, both procedures do not provide a good 
estimation on flapping acceleration. There appears to be a 
phase lag between the true and the estimated acceleration. 
The IBC Kinematic Estimator method utilises two 
accelerometers attached to different locations on the blade 
to extract both flapping angle and angular acceleration 
information. Algorithm is then used to estimate flapping 
angular rate via a kinematic observer. This technique 
provides good estimations of all three states used for IBC 
control law synthesis. The method assumes that the blade 
motion is due to rigid body flapping. 

The control for IBC can be exercised in the rotating 
frame by directly attaching actuators to the rotor blades. 
This mechanism, though appealing because of its simplicity, 
place severe demands on the actuators. It is shown that IBC 
can be achieved through a swash plate in the non-rotating 
frame when the control degrees-of-f reedom equals to the 
number of blades. This is the case for a three bladed 
helicopter. When the number of blades is four, coordinate 
transformation can still be utilized to find an equivalent 
swash plate feedback control. A simulation study with GENHEL 
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shows blade root pitch motion generated by IBC via swash 
pl.t, motion is reasonably close to that produced by direct 

IBC in the rotating frame. 

Also presented in the thesis is the general rotor 
states feedback control system in the non-rotating reference 
frame. Ham's IBC system is shown to be a special case of the 
general rotor states feedback scheme. The optimum state 
feedback gains in Equation 7.11 can be found for higher gust 
alleviation. This is beyond the scope of this thesis, and is 
not further pursued here. 
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APPENDIX A 


Derivation of Ham's IBC method 

The flapping equation of motion is described by : 

(3 + A ( U » 0 ) 8 + B(p,0)(3 = C(n,0)0 + Wg( U , W , Q ) (1) 

For the UH-60 blade used in the simulation, the parameters 
are given as : 

A = 23.76 ♦ 31 * u * sin0 

B = 734 + (692.24 + 1323.8 * ji * sin0) * u * cos0 

C = 684.3 + (1808 + 1313 * u * sin0 ) * 4 * sin0 

Wg = Kg (972 * sin wt + 792 * u (cos(Q-w)t - cos(Q+w)t)) 

e swp = 0.2975 + 0.009 * cos Qt - 0.142 sin Qt 

U =0.18 

Q =24 rad/second 

w =13 rad/second 

Here, the total blade root-pitch angle is consisted of two 
parts : 


6 



®ibc 
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where 8 SW p is the effective blade pitch input due to motion 
of the swash plate, and 0^ c is the individual pitch feedback 
control signal for gust alleviation. 

The feedback control law is generated as : 


e ibc = (KA (3/Q 2 + KR (3/Q + KP |3 ) 


( 2 ) 


where R’s are feedback control gains and 


® swp = * ®swp 


(3) 


K SW p is the forward controller gain. Substituting Equations 
(2) and (3) into Equation (1), we obtain 


A + C KR/Q B + C KP 

e + ( --) b + ( 0 = 


1 + C KA/Q 2 


1 + C KA/Q 


2 


) Kj W p 8 S yp + ( -y) Wg 

1 + C KA/Q 2 1 + C KA/Q 2 


(4) 


To reduce the gust effect acting upon the rotor blade, the 
controller gains can be selected such that the closed-loop 
blade dynamics equation will assume the form 


(3 ♦ A(u,0)(i ♦ B(U,*)P = C(u,0)0 swp + 1/KG H ff 
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The controller gains are assigned as : 


KR = KA * A/Q 

KP = KA * B/Q 2 

K swp = 1 + C * KA/Q 2 

KG = 1 + C * KA/Q 2 

For Ham’s simplified IBC scheme, the following assumption 
are made : 

A * Q ; B Z Q 2 ; C * Q 2 

With the above simplification, the controller gains become 

KR = KP = KA 
K swp = KG = 1 + KA 

The IBC feedback control law is now 
e lbc = ( (3 /Q 2 + B /Q + 3 ) k'A 
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APPENDIX B 


Derivation of Model Reference IBC method 

To design a IBC control law which is less sensitive to the 
DC value of the blade flapping motion, a Model Reference IBC 
system is studied. Flapping equation of motion is given as : 

S ♦ A(p,0)(3 + B(y,0)(3 = C(u,0)6 + W g (y,w,Q) (1) 

where 

0 = 0 s«p + 0 ibc < 2 > 

If a math model of the blade is known, 

(3 m + A(u,0)|3 m + B(u,0)3 m = C(w,0)0 swp (3) 

then the error in flapping angle can be formed by 
e ( t ) = (3 ( t ) - 0 m ( t ) 

D(t) is the measured angle, and P B (t) is the model response. 
Prom Equations (1) and (3), we have the error dynamics 
equation : 

e + A(y,0)e + B(U/ 0)e = C(u,0)e ibc ♦ H g (u,w,Q) (4) 
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A control law can be generated as 


0 ibc = (KA e/Q 2 + KR e/Q + KP e ) (5) 

where 

KR = KA * A/2 

KP = KA * b/Q 2 

Substituting KR and KP, Equation (5) is now given as 

0 ibc = KA/2 2 ( e + A (p,0)e + B(u,0)e ) (6) 

Equation (4) becomes 

1 

e + A(p,0)e + B ( p , 0 ) e = W. (7) 

1 + KA C/Q z 


Equation (1) is now 


13 + A(u,0)(3 ♦ B(u,0)|3 = C(u,0)6 sirp + W g (8) 

1 + KA C/2 Z 


Prom Equation (8) and Equation (4) of Appendix A, both 
control schemes reduce the effective gust level by a factor 
of 1 / ( 1 + KA ( C/Q 2 ) ). MRIBC is dependent only on the 
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1 - 1 '1 it is less sensitive to the DC 

error signal e, e and e. n is 

value of the flapping motion. 

From a point of view of practical application, it is 
desirable to have constant coefficients in the blade 
flapping math model. Using the same approximations as in 

Appendix A, 


A * Q 


B * Q 2 


C * Q 


2 


Equation ( 3 ) becomes 


+ Q + a 2 P* * ° 2 6 swp 


( 9 ) 


and the control law is 


’ibc 


= KA ( e/Q 2 + e/Q + e ) 


( 10 ) 
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Open loop 

IBC in effect 

Figure 4.1 Oust alleviation vs controller gain KA 
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Figure 4.3 IBC performance, Ham’s with time-varying gains 
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Figure 4.4 IBC performance. Ham's with time-invariant gains 
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Figure 4.5 Control input for small gain KA 
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Figure 4.8 IBC performance, MRIBC with time-varying gains 
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Figure 4.9 IBC performance, MR IBC with time-invariant gains 
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Figure 5.2 Direct Flap Estimator, with plant dynamics 
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Figure 5.3 Direct Flap Estimator, no plant dynamics 
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Figure 5.6 3 state Direct Flap Estimator 
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